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On The Spectrum of Yang-Mills Theory in 
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Abstract: We review our recent work on the glueball spectrum of pure Yang-Mills theory 
in 2+1 dimensions. The calculations make use of Karabali-Nair corner variables in the 
Hamiltonian formalism, and involve a determination of the leading form of the ground-state 
wavefunctional. 
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The understanding of the non-perturbative dynamics of Yang-Mills theory is one of the grand prob- 
lems of theoretical physics. The 2 + 1 -dimensional theory is expected on many grounds to share the 
essential features of its 3 + 1 dimensional cousin, such as asymptotic freedom and confinement, yet 
is distinguished by the existence of a dimensionful coupling constant. This simple fact has important 
consequences. 

Here, we describe the determination of the ground state wave-functional in a specific approxima- 
tion, the knowledge of which enables us to determine the mass gap, string tension and the glueball mass 
spectrum. The results are in excellent agreement with the lattice data in the planar limit [1]. Further 
details of this work, as well as a full set of references, may be found in [2, 3], Here we will concentrate 
on the main ideas and outstanding issues. 

We will describe the pure SU (N) Yang-Mills theory by transforming to corner variables [4]. One 
advantage of these variables is that the passage to a lattice regulator is straightforward. This change of 
variables may be done exactly, the Jacobian of the path integration measure being computable. Given 
a local coordinatization of space, we introduce the straight Wilson lines 



Mi (x) = Pexp 



(i) 



or equivalently Ai — —diMi M~ 1 ,\/i. Gauge transformations act linearly Mi i— ► gMi and local 
gauge invariant variables may be constructed, Hij(x) — M^(x)Mj(x). These are the (generally 
constrained) corner variables. In two spatial dimensions, a complex spatial coordinatization may be 
used, leading to the Karabali-Kim-Nair parameterization [5, 6] 



.4 



-dM M~\ A = M^ 1 dM\ H = M^M 



(2) 



These variables possess a new local holomorphic invariance 

M(z, z) ^ M(z, z)h)(z), M*(z, z) ^ h(z)M*(z, z), H(z, z) ^ h(z)H(z, z)h\z) (3) 
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which must be preserved. A holomorphic connection for this symmetry is 

J = dHH- 1 , J^hJh- 1 +dhh- 1 (4) 

and the corresponding covariant derivative D — 8 — J = M^VM* \ with V = 8 + A the usual 
gauge-covariant derivative. The path integral measure may be written 




The Hamiltonian, written as a functional differential operator in J's, was constructed by Karabali-Nair 
and takes the collective field form 



where m is the 't Hooft coupling m = g 2 N/2ir. In principle, wavefunctionals, expressed as functionals 
of J, may be found by solving the functional Schrodinger equation. Note that at weak coupling, the 
Hamiltonian is dominated by the trailing potential term, while at strong coupling, it is dominated by 
the kinetic operator. By looking at the action of parity and charge conjugation, it becomes clear that 
holomorphic invariant wavefunctionals may be constructed from 8 J — [D,8] and A = {D, 8}, 
= ^{dJ/m 2 , L — A /2m 2 ]. As we shall see, in considering the ground state wavefunctional, it 
appears to be useful to restrict our attention to the gauge and holomorphic invariant form 



We note that although this is not the most general functional (terms quartic in 8 J, etc., are allowed 
by symmetries), it is highly non-trivial in the sense that the kernel K(L) is allowed to be an arbitrary 
function, whose exact form will be determined by the Schrodinger equation, and its asymptotic features 
must be consistent with asymptotic freedom and confinement. We are assuming large N here, which 
implies a single trace but does not further simplify the exponent. The real importance of large N is 
in the expectation that interactions amongst gauge invariant states are eliminated and that the states 
are stable. Note further that the form (7) is not equivalent to a simple (covariant) derivative expansion, 
as all orders in derivatives are included. Since the magnetic field B = M^dJ , it appears that 

this may be thought of as a sort of gauge curvature expansion. In a sense, 8 J may be thought of as 
an adjoint constituent, with glueball states, which are normally identified with vibrating closed strings, 
modeled as configurations of a pair of constituents connected by strings. Whether or not this picture is 
borne out, is a matter for experiment to decide. As we detail below, we are able to compute glueball 
masses which agree remarkably well with the available lattice data. The form (7) may not be sufficient 
for the calculation of other observables. 

The first term in the kinetic operator counts the number of J's in a functional; through detailed 
computations, it appears that the second term in T acts essentially to restore holomorphic invariance. 
This implies that, defining O n = tr 8JL n 8j, 



The ellipsis contains mixing with operators containing more factors of 8 J which will not concern us. 
Eq. (8), crucial to the results that follow, may be explicitly demonstrated for low values of n, but has not 
been firmly established in general. The required calculations rely on a holomorphic invariant regulator 
and the calculations are tedious; better calculational methods remain to be found. 

Now, given the assumed form of the ground-state wavefunctional (7), the Hamiltonian (6) and the 
result (8), we find that the Schrodinger equation takes the form 



Hkn[J] =T + V = m 




(7) 



TO n ~ (2 + n)O n + ... 



(8) 




(9) 
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The first ellipsis contains (divergent) terms contributing to the vacuum energy, while the trailing ellipsis 
contains terms of higher order in 8 J. Given the assumed form of the vacuum wavefunctional, this 
truncation of the Schrodinger equation is appropriate and consistent. The quantity 1Z may be found by 
regarding K as a power series in L, and one finds 

K = -K(L)-^K(L)+LK 2 (L) + 1 (10) 

The Schrodinger equation requires that we set this to zero, resulting in a differential equation for K 
of the Riccati type. Through a series of redefinitions, this may be cast as a Bessel equation, and one 
obtains 

K { L)=**&®. (11) 

Although there are other solutions to the differential equation, no other is normalizable in the given 
path integral measure. We note that although this is a very complicated function (which we now regard 
as a function of momentum), it has the asymptotics K(p) — > 2m /p as p — > oo (asymptotic freedom) 
and K{p) — > 1 as p — > (confinement). We have thus found that the only normalizable ground state 
wavefunctional is consistent with both confinement and asymptotic freedom. As argued by Karabali 
and Nair, the infrared limit of our wavefunctional K implies a string tension yfo ~ ^= obtained 
by a dimensional reduction argument. This result agrees precisely with lattice data. However, it is 
not possible by this line of reasoning to detect different string tensions appropriate to different gauge 
representations. It is believed that this is consistent, being an artifact of the large N approximation 
known as Casimir scaling [7] (in particular, this result would be inconsistent at finite N). In effect, 
there is an order of limits problem at large N and large distance. 

The ratio of Bessel functions has a rich analytic structure and encodes the mass spectrum of the 
theory. By Fourier transforming, we find 

oo 

K~\\x y\) = - , J T{M n )We- M «\*-v\ (12) 

^2n\x - y\ ^ 

where M n — 7 2 ,„m/2 and J2(72,n) = 0. As we will see, this result has direct consequences for the 
mass spectrum. 

To probe the mass spectrum, we consider pair correlation functions of gauge invariant operators 
with definite spin, parity and charge conjugation quantum numbers. The simplest ++ probe operator, 
for example, is trdJdJ, and we wish to compute 

(trdJdJ(x)trdJdJ(y)) = J dfj,[H]e 2NSwzw[H] |* | 2 trdJdJ{x)trdJdJ{y) (13) 

at large spatial separation. To proceed, we first rewrite the measure as an integral over J. Note that if 
we introduce a variable J = dHH -1 , then we have the 'reality condition' 8 J — [D,J]. Given the 
gauge transformations 8 A ~ [V, 5MM~\5A ~ [V, M^8M^\ , the measure is 

with dfi[H] the measure corresponding to dsf nv — J Tr(8HH~ 1 ) 2 . 

We would like to now transform this to the J variables, namely, we would like to find the measure 
corresponding to the distance ds 2 j = J Tr(SdJ) 2 . We first note that 8dJ ~ [8, [D, SHH^ 1 }}. So we 
have 

dfj,[A] _ detVV d^i[3j] 
Vol G ~ DetdD Vol G ho i 
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Since D = M'VM"' and 3 — VM~^, the determinant factor formally cancels exactly. It should 
be appreciated here though that this is a formal result: in particular, although not apparent, this is a 
finite volume measure, because the integration is performed over a finite volume region in the complex 
8 J space. Ignoring such subtleties we obtain 

(trdJdJ(x)trdJdJ(y)) = J dfi[dJ}\^ \ 2 trdJdJ{x)trdJdJ{y) (15) 

To the approximation in which we ignore interactions (that is, take K(L) to be a function of momen- 
tum), we may thus regard d J as a 'constituent,' with correlators determined by K(p). Consequently, 
we find 

(trdJdJ(x)trdJdJ{y)) ~ R- 2 (\x - y\) = ^ -#_ e -{M n +M m )\ x - y \ (16) 

m,n ' X Ml 

a form consistent with a single particle ++ pole of mass m m ,n = M m ± M n = (72,™ ± 72, n), m/2. 
We note that K^ 1 is not the propagator of a physical mode, but it does determine physical propagators. 
The resulting masses for the ++ are collected in the table and compared to the lattice data [1]. The 
numbers quoted are in units of the string tension. 



State 


Lattice, N — > oo 


Our prediction 


Diff, % 


0++ 


4.065 ± 0.055 


4.098 


0.8 


0++* 


6.18 ±0.13 


5.407 




0++" 


6.18 ±0.13 


6.716 




0++*** 


7.99 ±0.22 


7.994 


0.05 




9.44 ± 0.38 


9.214 


2.4 



We see very good agreement, apart from the second and third states. We note though that the 
average of these two states coincides with the second lattice state, and we may take this as a prediction 
that these two states were not resolved in the lattice studies. 

Similar results may be obtained for other states using suitable probe operators and results are 
equally encouraging. The states for example are given in the second table. 



State 


Lattice, N — > oo 


Our prediction 


Diff,% 


0"" 


5.91 ±0.25 


6.15 


4 


— , 


7.63 ± 0.37 


7.46 


2.3 


** 


8.96 ±0.65 


8.73 


2.5 



The resulting spectrum appears to have an exponentially rising density of states, which may be 
taken as a manifestation of the effective QCD string. 
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